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i . We point out that for a large class of parametrized theories, there is a constant in 
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the constrained Hamiltonian which drops out of the classical equations of motion in 

configuration space. Examples include the mass of a relativistic particle in free fall, 

the tension of the Nambu string, and Newton's constant for the case of pure grav- 
en ' 

ity uncoupled to matter or other fields. In the general case, the classically irrelevant 
constant is proportional to the ratio of the kinetic and potential terms in the Hamilto- 
nian. It is shown that this ratio can be reinterpreted as an unconstrained Hamiltonian, 
which generates the usual classical equations of motion. At the quantum level, this 
immediately suggests a resolution of the "problem of time" in quantum gravity. We 
then make contact with a recently proposed transfer matrix formulation of quantum 
gravity and discuss the semiclassical limit. In this formulation, it is argued that a 
physical state can obey a (generalized) Poincare algebra of constraints, and still be an 
approximate eigenstate of 3-geometry. Solutions of the quantum evolution equations 
for certain minisuperspace examples are presented. An implication of our proposal is 
the existence of a small, inherent uncertainty in the phenomenological value of Planck's 
constant. 
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1 Introduction 



Actions which are invariant under reparametrizations of the time variable typically have 
Hamiltonian constraints of the form H[p,q] = 0, which can be viewed, in an initial value 
problem, as constraining the initial state {p, q}^. In the Dirac quantization scheme, 
these constraints on the conjugate variables become constraints on the physical states, 
so that instead of the usual Schrodinger evolution equation one has 



H 



-i— Q° 



q[q a ] = 



In cases where the constraints are parabolic, such as the parametrized non-relativistic 
particle, or parametrized scalar field theory, the constraint equation (TO) can often be 
treated exactly as a Schrodinger equation; the coordinate whose derivatives appear only 
to first order is identified with time (in the non-relativistic case), or "many-fingered 
time" (in the scalar field-theory case) [TJ. On the other hand, for hyperbolic constraints 
such as the Klein-Gordon equation with an arbitrary background metric g^ u 

d 2 
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or the Wheeler-de Witt equation 
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it has proven difficult to identify an appropriate evolution parameter and a unique, 
positive, and conserved probability measure. In quantum gravity, this difficulty is 
known as the "problem of time"; c.f. ref. || and 0, for recent reviews. 

In this article we propose an alternative to the conventional procedures for quan- 
tizing parametrized theories. Our proposal begins with the rather trivial observation 
that, at the classical level, there is no observable difference between an action S and the 
same action multiplied by a constant. Since, e.g., the action for a relativistic particle, 
the Nambu string, and the Einstein-Hilbert action are, respectively 
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dr dr 
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it follows that the mass m, the string tension T, and Newton's constant Gn = k 2 /167t, 
which appear in the Hamiltonian constraints of the relativistic particle, string, and pure 
gravity theories respectively, do not appear in the Euler-Lagrange equations, which are 
the geodesic equation 

for the relativistic particle, the equation of motion 

V ^ X ^W) = (6) 
for the Nambu string (in orthonormal coordinates), and Einstein's equations 

= o (7) 

for pure gravity. The mass of a particle (or the tension of a string) can never be 
determined from its trajectory in free fall, neither can Newton's constant be determined 
from vacuum solutions to Einstein's equations. In the theories described by (f|) these 
are classically indeterminate parameters, which only become relevant by introducing 
non-gravitational fields and forces (i.e. changing the theories). 

In section 2, below, we generalize this observation to any time-parametrized theory 
with hyperbolic constraints, in particular, to gravity coupled to other fields. It will be 
shown that there is always a classically undetermined parameter £, which is equal to 
the ratio (denoted JE) of the kinetic and potential terms in the constrained Hamiltonian 

c mi i Kinetic[p,g] 
- g = ^g] = PotentiaI[ g ] (§) 
We will also show that replacing the usual (constrained) Hamiltonian H[p, q] in Hamil- 
ton's equations by the expression JE\p, q] generates, in configuration space (superspace), 
the usual classical solutions of the theory. 

Our proposal for resolving the time problem in quantum gravity is based on this 
reasoning: since £ is indeterminate at the classical level, we see no reason that it must 
be regarded as a fixed (bare) parameter at the quantum level. Rather, let — £ denote 
the possible eigenvalues of the operator JE. Because JE ^ and, as shown in section 2, 
the Poisson brackets 

d T Q = {Q,JE} (9) 

generate dynamics at the classical level, it is proposed that JE is also the appropriate 
evolution operator to use at the quantum level, i.e. 

d T < Q >= - < Uq,JE] > (10) 



The Hilbert space of physical states is spanned by stationary states 



(11) 



and this time-independent Schrddinger equation will be recognized as just the usual 
Hamiltonian constraint with a given value of the parameter £ , which varies among 
stationary states. 

In section 3 we rederive the JE evolution operator in a completely different way, 
based on a transfer-matrix quantization of parametrized theories. This section is a 
review and extension of the approach suggested recently by one of us in ref. j3j, 
which leads to the same picture as that advocated in section 2. The transfer matrix 
approach has the advantage of fixing the operator-ordering in JE[p, q] and the functional 
integration measure, for parametrized theories with a discrete number of degrees of 
freedom. 

Section 4 is devoted to a discussion of the semiclassical correspondence between 
(H) and ([TUP, and in particular the relation between the evolution parameter r that 
appears in those equations, and the "many- fingered time" of general relativity. We 
argue that the space of physical states includes states which are sharply peaked (at a 
particular value of r) around a given 3-geometry. In such states, the dispersion around 
the given 3-geometry is, crudely speaking, inversely proportional to the dispersion in 



We also argue that a dispersion in S, which is characteristic of non-stationary states, 
would appear experimentally as an inherent uncertainty in the phenemonological value 
of Planck's constant. Unfortunately, we are unable to place on this uncertainty a 
reliable lower bound. 

In section 5 we illustrate our formalism by solving a simple minisuperspace model, 
and display the "wavepacket of the Universe" evolving through the stages of expansion 
of the scale factor, recollapse, and reexpansion. Section 6 contains some concluding 
remarks. 



We consider first a time-parametrized theory with D degrees of freedom, described by 
the action 



2 The Kinetic/Potential Ratio 
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G ab p aPb + m V(q) 



(12) 
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where the " supermetric" G a b has Lorentzian signature { — h +...+}. For V(q) =const., 
this action characterizes the motion of a relativistic particle in the spacetime metric 
G a b, other choices of V arise in minisuperspace models of quantum gravity. Making 
the rescaling p a — > p a /\fS (an extended canonical transformation) the action becomes 



= —=S[m = VSm ] (13) 

It is clear from this equation that the trajectory q a (t) in configuration space connecting 
two arbitrary points q$ and cannot depend on the value of m , which is therefore 
classically indeterminate. We may also write 

1 



where 



V£ S£ (14) 



S £ = Jdt( Pa ^-NH £ ) 



H e = 2v ^ m G ab PaPb + VSm V (15) 

Suppose S is stationary for a trajectory in phase space q a {t) = Q a (t), p a (t) = 
P a {t), N(t) = A/"(t), which connects the generalized coordinates q$ and qf. It fol- 
lows that S £ is stationary for q a {t) = Q a (t), p a {t) = VSP a (t), N(t) = Af(t). This 
means that the trajectory in configuration space q a (t) = Q a (t) connecting and q^ 
will be the same whether the fundamental action is taken to be S or S . It is in this 
sense that the parameter S is indeterminate, since trajectories in configuration space 
generated by 

dt dp a dt dq a y 1 

will be independent of E. 

Let us rewrite the H £ = constraint as 

-£ = M[p,q}= ^ G l aVh (17) 

It is easy to show that the configuration space trajectories generated by treating JE[p, q] 
as though it were a Hamiltonian 

M _ g dq a _dM dp a dM 
dr dp a dr dq a 



are equivalent to the solutions generated by The first of these equations, M 

—£, is just the Hamiltonian constraint 7i £ = 0. The second equation is 



dq a 
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and for the third equation we have 
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dp a 

m V dp a 
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dq a 2mo 



T-G ab p a p h d 



moV dq a 







m V dq a 



G ah p a p b + £m V 
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01), and M 



m V dq a 

Suppose q(r), P (r) is some particular solution of ([19] 
choosing 

NiT) ^ vm 

and rescaling r by \f£/rriQ to give r the conventional units of time, we see that this 
is also a solution of the Hamilton's equations and constraint ([IB]). It follows that, in 
general, the Poisson bracket evolution equation 



(20) 
-£. Then, 

(21) 



d T = {O, M} 



(22) 



supplemented by JE = —£ will generate time evolution which is equivalent, up to a 
time reparametrization, to evolution generated by 



d t O = {O, NH £ } 



(23) 



supplemented by TC £ = 0. 

The extension of these remarks to general relativity is fairly straightforward. Let 
us denote the action of a generally covariant field theory as 

dq a 



S = J d^x[p a ^-N7i x -N t 7i} x _ 
H x = K 2 G ab p a p h + ^U{q) 



(24) 



3 It should be noted that, since M is dimensionless, the evolution parameter r in eq. (18) has units 
of action. 



where the q a (x) are the set of all fields, gravitational and non-gravitational, G a b is the 
supermetric, N and iVj the lapse and shift functions, and TC X are the supermomenta, 
linear in the canonical momenta p a {x). It is convenient to rescale all non-gravitational 
fields by an appropriate power of k so that all q a (x), and all G ab , are dimensionless. In 
the special case of pure gravity, we would identify 

{a = 1 - 6} <-> i < j} 

q a (x) «-> g i:j (x) 

P a (x) «-> 

Gabix) 



p 13 (x) 
2p ij (x) 

x) 



(i < J) 



£<ijnm 



u 



1 



- 3 R) 



(25) 



By the same rescaling p a 
level to the action 



P, 



,/ y/S as before, the action S is equivalent at the classical 



d 4 x[Pa^ ~ NH £ X - N t Hl 



H £ x = ^=G ab PaPb + V£^U(q) 
Variation of S £ with respect to the lapse gives the Hamiltonian constraint 



(26) 



Hi = ^G ab PaPb + Vs^gU(q) = (27) 

which, multiplying again by an arbitrary function N(x), we write in the form 

Jd 3 x NK 2 G ab PaPb 

Jd*xN^gU(q) [ } 

valid for any N(x). Because of the supermomenta constraints 7i x = 0, we can add a 
term proportional to the supermomenta to the ratio in fl28|) without changing its value. 
Define 



M\p,q,N,Ni\ 



d x 



Nn 2 G ab p a p b 1 



/ d 3 x' y/gNU (q) m P 

— I d 3 x {N^G^PaPh + NiHi} 
ip J 1 > 



nip 



where 



N(x) 



mpN(x) 



Jd^x' ^gNU(q) 



(29) 



(30) 
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and mp is an arbitrary mass parameter. We will now show that M can itself be rein- 
terpreted as a Hamiltonian, in the sense that the classical orbits in configuration space 
obtained from the corresponding Hamilton's equations are identical to the extrema of 
the action (|24). 

The equations of motion derived from M are 



dq a (x) 
dr 



dpa(x) 
dr 



5M 

— f d 3 x' 
mp J 

5M 



N(x')-^—(K 2 G cd PcPd ) x , + N t (x')^-V\, 



5pa(x) 



5p a (x) X 



X 



d 3 x 



mp ^ 
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fN^gU 
d 3 x' 



d 3 x' 
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mpN ^Sr{xj^ U)x ' 
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mp 
+N t (x' 

5M 
6N(x) 

1 

JNy/gU 
5M 1 



-7~C 

Sq a (x) x ' 



(K 2 G ab PaPb ) a 



(jNK 2 G ab PaPb 



(31) 



5Ni(x) mp 

The value of M is of course a constant of this motion, denoted JE = —£, and r is 
dimensionless. From the definition of in eq. (|26|), and rescaling t = y/Sr/mp (to 
give t the dimensions of length), the equations of motion ([H]) become 



dq a (x) 

dt 
dpg(x) 

dt 



d 3 x' 
d 3 x' 



N{x')^r^n% + Ndx')^-U\, 



5p a (x) 



5p a (x) 



n_/i 



K 

7s 





5q a (x 

G ab PaPb + VS^U = 



5q a (x) 



(32) 

where the shift function is Ni = Ni/yfS. Comparison of the equations of motion (R5J) to 
the Hamiltonian equations of motion that would be derived from S of eq. fl26|) shows 
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that they are the same equations, apart from the restriction of the lapse function iV in 
(p^) to satisfy 

J d 3 x ^fgNU = m P (33) 

The trajectory in configuration space is independent of the choice of lapse; it is also 
independent, as noted above, of the choice of 8. 

We can therefore conclude that JE\p, q, N, iVj, viewed as a Hamiltonian, generates 
the same classical dynamics as the more conventional Hamiltonian form 

H = J d 3 x[NH x + NiTQ (34) 

The crucial difference between M and H is that H is constrained to vanish, in the 
standard formulation, whereas the value of M is unconstrained: the constant £ can take 
on any value on an orbit, depending only on the initial choice of {p, q}. Obviously, this 
difference is quite important at the quantum level. Instead of the usual Wheeler-de 
Witt constraint equation 

HV[q] = (35) 

we will obtain a Schrodinger equation 

MV[q,T]=ihd T V[q,T] (36) 
Several questions arise immediately: 

1. If (^) is the dynamical equation of quantum gravity, what has become of the 

Poincare algebra of constraints which represents the diffeomorphism invariance 
of the theory? 

2. What is the relationship of the Schrodinger equation (|3~6"D, based on the iE-operator, 

to path-integral quantization? 

3. What is the operator ordering in M\p, q, N, iVj], and what integration measure should 

be used for the inner product of states? 

4. The time-evolution parameter r in the Schrodinger equation is only a single variable; 

how is it related to the "many-fingered time" of general relativity? 

To answer the first question (the others will be dealt with in subsequent sections), 
let us note that physical states ^/[q, r] in the Schrodinger representation must be inde- 
pendent of the functions N and iVj. Expanding an arbitrary \1/ in stationary states 

nq,r]=J2ae^[q]e t£r/h (37) 
e 



o 



where 

m$ e [q] = -£$e[q] (38) 
the condition of Aj-independence gives 

5 

oNi(x) 

= K*e (39) 

while ^-independence requires 



= M®** 6 



' r K 2 G ab PaPb - y/gUM} x $ e 



Jtfx'^/gNU 



H\$e (40) 



where we have provisionally taken an operator-ordering in M with momenta to the 
right, and also used the supermomentum constraint (|55|). Therefore, the Hilbert space 
of physical states consists of linear combinations of iE-eigenstates 

( 41 ) 

£ 

each of which satisfies the constraints 

H e x ®£ = -^=[K 2 G ab p aPb + £^/gU] x $ s = 

Hi$e = (42) 

Apart from the parameter £ , these are the standard constraints on physical states 
of quantum gravity. As is customary in this subject, we assume that there exists 
some operator-ordering such that the commutators of operators H x , 7i x close on the 
Poincare algebra, as do the Poisson brackets of the corresponding classical quantities. 
Given this assumption, the physical states of the form (51) satisfy, in our formulation, 
the generalized Dirac constraints 

{[K 2 G ab p a p b } x -[^g-U] x m}y phys = 

K% hys = (43) 

In the case of pure quantum gravity, the 8 parameter has a simple interpretation. 
The Hamiltonian constraint becomes 

j^G^pV' + ^v^(- 3 ^)} ®£ = (44) 

i n 



and therefore 



H 2 



G " = Tw? (45) 

is the effective Newton's constant for the degenerate subspace of physical states satis- 
fying iE$ = — The space of physical states is spanned by states $g which satisfy 
the standard Hamiltonian and supermomentum constraints, but with a different value 
of Newton's constant associated with each S. In general, physical states are not eigen- 
states of Newton's constant. 

To summarize our proposal for pure gravity: since Newton's constant is indetermi- 
nate at the classical level, we see no reason that it is necessarily a fixed bare parameter 
at the quantum level. Instead, treating Gat as a quantum number allows for a vast ex- 
tension of the space of physical states. It is this extension which we propose to exploit, 
as seen in eq. (PE|), to resolve the time problem in quantum gravity. In the general 
case of gravity coupled to other fields, it is not Newton's constant per se but rather the 
kinetic/potential ratio in the constrained Hamiltonian that is classically indeterminate. 
Regarding this ratio as a q-number leads again to an extension of the space of physical 
states and, as we will argue further below, to a solution of the problem of time. 



3 The Transfer- Matrix Formulation 

We now give an alternate derivation of the Schrodinger equation (|36|) for parametrized 
systems, following the transfer- matrix formulation proposed by one of us in ref. 0. 
The advantage of this approach, apart from making the connection to path-integral 
formalism, is that it also fixes the measure and operator-ordering, at least for systems 
with a discrete number of degrees of freedom. 

The transfer matrix T e for non-parametrized systems, in statistical mechanics and 
Euclidean quantum mechanics, is an operator which evolves states by a time-step 
At = e 

= Jd D qi x(q)ex P [-S t (q',q)}^(q,t) (46) 

Denoting by S[q2,qi] At] the action of a classical solution q(t) running between the 
initial point q\ at time t and final q2 at time t + At, the expression S e in Euclidean 
quantum mechanics is given by the continuation of 5* to imaginary time lapse 

Se(q2, qi) = iS[q 2 , qi, ie]/?i (47) 
The measure /i(q) is defined such that T e is an identity operator as e — > 0. The usual 



i i 



quantum-mechanical Hamiltonian is essentially the logarithm of the transfer matrix 

H = jim(--) H%] (48) 

and the continuum (Euclidean) path integral is defined as the limit of a product of 
transfer matrices 

V(<Z', *i) = / Dq(t < t < tJe-^V^iq, t Q ) 

.N-1 N-1 
= l }j$J II d °Qn /i(g n )exp[- J2 S t (q n+1 ,q n )}ip(q ,t ) 



lim(T e ) N M,t ) where N = - (49) 



Now let us again consider parametrized theories of the form 

dq a 



S = Jdt( Pa -^-NH) 



U = J-G ab PaPb + mV(q) (50) 
2m 

For parametrized theories of this sort, the transfer matrix formalism above breaks 
down at eq. (f47l). The problem is that the action S[q2, qi, G a b] of a classical trajec- 
tory connecting initial coordinates q\ and final coordinates q 2 depends only on those 
coordinates (and choice of supermetric G a b); there is no additional dependence on a 
time lapse At, and the values of parameters T\ and T2 that happen to be associated 
with the initial and final coordinates are irrelevant. A transfer matrix 7^ based on P7| ) 
would then be independent of e; this is one of the obstacles encountered in trying to 
apply directly the path-integral approach to find the time-evolution of physical states 
in parametrized theories. 

It is premature to conclude, however, from the failure of (|47|) , that the transfer 
matrix concept is meaningless for parametrized theories. Let us denote an orbit in 
Hilbert space by ip(q,T), and require that the transfer matrix T t evolves states in the 
parameter r such that the orbit of the center of the wave packet 

<q a (T)>=<iP(q,T)\q a \iP(q,T)> (51) 

obeys an appropriate Ehrenfest principle. We will show, following ref. that this is 
achieved by replacing (|47|) in the definition of the transfer matrix by 

S e (< b ,q 1 ) = iS[< b ,q 1 ,G% > ]/Vdi (52) 

where the Euclidean rotation At — > iAt in (f47D is replaced by a rotation of the signature 
of the supermetric 

G ab (q) = El{q) Vl3 Ei{q) 

G E ab = ^[V{q)]E^q)5 l3 El{q) (53) 
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The reason for this peculiar rotation of signature is to ensure that S e (q2, qi) is real 
for all choices of q\ and q 2 . To show this, we first compute S[q, q + Aq, Gab], beginning 
from Hamilton's equation 

q" = N^ = -G ab p b =► p a = -G ab q b (54) 
op a m JM 

Insert this into the constraint equation TC = 

' G ab q a q b + V = (55) 



2N 2 

and solve for the lapse 



N = [~^G ab q a q b ] 1/2 (56) 



Then 

rAt 



S[q,q + Aq,G ab ] = / dt -G ab q a q b 

Jo iV 

= -m [ dt J-2VG ab q a q b 
Jo 



mJ-2VG ab Aq a Aq b 



= -m^-g ab Aq a Aq h (57) 

where we define a modified supermetric 

Q ab = 2VG ab (58) 

It is clear that if the signature of G ab is Lorentzian, then S[q, q + Aq, G ab ] is not 
necessarily real. However, upon the signature rotation (|53|), the function 



S e (q,q + Aq) = mJ2VG^ b Aq a Aq b /Ve?i 



mJggAq a Aq b /Veh (59) 



is strictly real- valued for any choice of Aq a . 
The integration measure fi is chosen to be 



H~\q!) = (V7h) D lim / -^3- eMSe(q', ?)) (60) 



eh) 

This choice ensures that T e — > 1 in the e — > limit, and that the symmetries of the 
action are reflected in the measure. 

The equations (|52D,(|53D,(|50T), when inserted into eq. (^§), define a transfer matrix 
T e for the parametrized, signature-rotated theory, which in general depends on Q E and 



1 o. 



det [Q E ] . The final step is to undo the signature rotation in computing the corresponding 
evolution operator M 



hrn(-~)ln[T e ] 



g E ^g 

dct(gs)^| det (g)| 



(61) 



and the evolution of states is given by 

ihd T tp(q, t) = Mip(q, r) 



(62) 



We can now evaluate T t and M for actions of the form ([5(]), leading to the function 
S e shown in (^). Begin by introducing Riemann normal coordinates £ a around q' a = 
q a — Aq a , which bring Q® h = 5 a b at £ a = 0. In these coordinates 



S e (0, = mJ5 ab ^e + 0(e) /Veh 



(63) 



The 0(£ 5 ) terms will not contribute to M in the e — > limit, and can be dropped. The 
measure is then 

, r ,j D p BAn^ , 

Vr\4) = (v^Jimy jjjLj det[-^-]«p[-m|e|/V^] 



i 



a 



T(D/2) ' ' m ' jQEtf) jgEtf) 

The operation of the transfer matrix ( [46]) becomes 

d D Aq 



(64) 



^(q',r + e) = J -y/G E {q' + Aq) exp[-my/g^Aq a Aq^ /V^]ip(q' + Aq, r) 

^ (1 _i^ + ... )exp[ _^| ] W T) + ^ 



i + eh^-^d^d,, - th^X^n + ou 2 ) 



(65) 



2m 2 6m 2 

where TZ is the curvature scalar formed from the modified supermetric Q E b of eq. fl58D , 
which has been rotated to Euclidean signature. Transforming back from Riemann 
normal coordinates, we have 



% = !+€ 



n- 



2m 2 ^(p dq n 
from which we extract, according to (|6| 



dq Ti 



Qm 2 



+ 0(e 2 



D + l. 2 1 d 
-w 



2m 2 



Q I dq n 



dq r 



Qm 2 



(66) 



(67) 



i i 



where Q = det(£? mn ), and G mn has been rotated back to Lorentzian signature. 

The classical quantity JE d corresponding to the operator JE is obtained by replacing 
derivatives with c-number momenta 

JE cl [q a ,p a ] = hm Mg a , - Pa ] (68) 

h^o oq 

which gives 

JE d = (D + 1) 2m l aVh (69) 
mv 

This quantity, apart from an unimportant overall factor of (D + 1), is simply the Ki- 
netic/Potential ratio introduced in the previous section. The transfer matrix approach 
is thus a second way of obtaining the evolution operator JE, and also provides a definite 
prescription for the integration measure (fi(q) oc \det(Q)\), and operator ordering. 

As a further consistency check, we note that in ordinary Euclidean quantum me- 
chanics we may define 

pip(q') = lim / d D qjj,m— — exp[— SJq', q)]ip(q) 



S e (q',q) 



1 {q'-qf 

-m h V(q)e 

— C 



/h (70) 



from which it is easy to evaluate the commutator 

[q,p]ip = lim J d D q fi{q'm— — — m— —q} exp[— S e (q, q)]ip{q) 

= hip (71) 

The factor Ti, rather than ih, is due to the Euclidean time. To derive the corresponding 
result in the parametrized theory, we begin from Hamilton's equation 

dq a dJE d 
dr dp a 

+ =* » - ^ C (72) 

- ^r G Pb Pa ~ DTi Gab ^ (72) 

and define 

r m 2 V <n' b — n b ) 

pM) = limy d D q Mg)p^T)^ - q ) ex V [-S e (q',qmq) (73) 

with S e and //(g) given by ( fa^f ) and (pCf). It is then easy to show that 

[q a ,p b ]^ = h5^ (74) 

as in the non-parametrized case. 

1 K 



Next we consider general covariant field theories of the form (|24]) . The main obstacle 
to computing S[q', q, Gab], and hence to computing the transfer matrix, is the presence 
of the shift functions Ni in the classical action. Covariant derivatives of the shift 
functions appear in Hamilton's equation, e.g. in the case of pure gravity 

dgij 



Of 



2K z NG ijnmP nm + N i;j + N r< 



(75) 



and these make it impossible to solve for the lapse-shift functions algebraically, in 
terms of dtQij- To determine the lapse-shift functions, it is necessary to solve certain 
intractable partial differential equations. Instead, we adopt the strategy of simply set- 
ting Ni = 0. In that case, of course, the supermomentum constraints are not obtained 
by extremizing the action, and must be recovered by imposing operator constraints of 
the form 

Q x [q a ,Pb]V =0 (76) 

on the physical states. These constraints (up to operator-ordering ambiguities) will be 
determined below. With the choice Ni = 0, Hamilton's equations give 

1 



Pa 2k*N 
Inserting this into the constraint equation 



G ab q b 



o = n x = j^G ab q a q b + Vou 

and solving for the lapse, gives 



(77) 
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N 



.1/2 



G ab q a q h 



(79) 



so we have 

AS = S[q',q,G ab ] 



- [d 3 x dt^-g^UG ab q a q b 

K J JO 

- f d 3 x[^-(g l >UG ab ) Aq»Aq» + 0(Aq 2 



dMV9)oW-(Ga b ) Aq"Aq b + 0{Aq 2 )] (80) 



where Aq a = q a — q' a , and we define 



Gab = ~~^U G ab 



(81) 



The notation (..)o means that the quantity in parenthesis is to be evaluated at Aq = 0; 
i.e. at q'. 



Applying the prescriptions (j46|) , (p2|) , (|53|) , the transfer matrix is formally obtained 
from 



if>tf, T + e) = J Dq fi(q)e~ AS/V ^ ifj{q') + J d 3 x 

5 2 if) 



6i/j 



5q a (x)8q b (y) 
^(g / ,r)+e[To + T 1 + T 2 ]+0( e 2 ) 



5q a (x 
Aq a (x)Aq b (y) + 



Aq a {x) 



(82) 



where the T n represent terms with n functional derivatives of if). To find these terms, 
we need to evaluate 



< Aq a (x 1 )Aq\x 2 ) > = J D(Aq) (fi) Aq a (x 1 )Aq b (x 2 ) 

1 



x exp 



Kj J 



(83) 



Unfortunately this quantity is highly singular, and in fact ill-defined without a regu- 
larization procedure of some kind. 

The authors are not aware of a non-perturbative regulator of the integral over 



three- metrics in (82) which preserves an exact diffeomorphism invariance. In the ab- 
sence of such a regulator (which is also crucial for sorting out the operator-ordering 
issues JT|), we can only make some general remarks about the regularized form of 
< Aq a (xi)Aq b (x 2 ) >■ Inspection of (|83D shows that < Aq a (xi)Aq b (x 2 ) >= for 
X\ 7^ x 2 ; we also see that (Q^ b )oAq a (x)Aq b (x) transforms like a scalar. One therefore 
expects the regulated expectation value to go like 



< Aq a (x l )Aq b {x 2 ) > 



reg r 



(Q Eab ^ 



V9 



) S 3 (x 1 - x 2 ) 



14) 



where v is a scalar quantity with dimensions of volume, v — > as the regulator is 
removed. A very important issue is whether v has some dependence on the 3-metric 
(gij)a, and this depends on the properties of the (unknown) regulator. If, as is the case 
with dynamical triangulation, there is a fixed short-distance cutoff I, then we expect 
that v ~ I 3 , and the number of degrees of freedom changes with the volume of the 
manifold. If, on the other hand, the number of degrees of freedom N p is fixed (as in 
the Regge lattice), then the volume per degree of freedom changes with the volume of 
the manifold. In the latter case, it is reasonable to expect that 



< Aq a (x 1 )Aq b {x 2 ) > 



reg r 



V 



1 

V9' 



15) 



1 7 



where V is the volume of the 3-manifold described by (gij)o, and (3 is proportional to 
the number of degrees of freedom in the regularization. A naive lattice regularization, 
replacing, e.g. 

Aq a (x) <-> Aq a {n) 
V 



"P n=l 

5 NpJ^Hj d 



5q a {x) V dq a (n) 

Dq «-> l[d D q(n) (86) 

n 

does, in fact, lead to (p5|) (c.f. [§J), although of course such a regularization does not 
at all respect diffeomorphism invariance. We will assume that there exists some good 
regularization leading to ( p5|) , although our justification for this assumption is largely 
a posteriori. 

Using eq. (BH), the T2 term (in (B2f)) contributing to the transfer matrix is easily 



evaluated; it is the only term which is important in computing the semiclassical limit. 
The other terms are operator-ordering contributions which, in the absence of an explicit 
regulator, we will ignore. The corresponding JE operator is 

»m» = -h 2 ^- [ d 3 x u- l G ab -^— (87) 

V J 5q a 5q b v ' 

where " JE" is in quotes to emphasize that this operator, by itself, does not yield the 
correct equations of motion in the classical limit. This was to be expected, because 
JE has been derived by setting the shift functions iVj = 0, a step which requires addi- 
tional constraints on the physical states. Now note that because of the spatial volume 
denominator in (0), we can write the evolution equation (|36"D as 



d 3 xQ x ^ = (88) 
where 

Q x = -tffU-iG*^ - ihjg± (89) 

and where we have absorbed the constant (3 in " JE" into a rescaling of the evolution 
parameter r in eq. (|36|). The extra constraints which need to be imposed on the 



physical states, which then generate the usual constraint algebra of general relativity, 
are simply 

Q x ^ = (90) 
at every point x. To show this, consider an arbitrary solution of the evolution equation 

^(q,r)=J2aee i£T/h Mq) (91) 
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Since the ag are arbitrary, the Q-constraint (pOf ) requires that for each stationary state 
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+ V£y/gU U g = 



(92) 
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But this is simply the Wheeler-de Witt equation for the action Sg of eq. (|26|)! Moreover, 
the constraints 7Yf $g = imply, via the Moncrief-Teitelboim interconnection theorem 
P], that the supermomentum constraints 



HI®, 







(93) 



are satisfied as well. In this way, we find that the stationary states {$£•} satisfy the 
usual constraint algebra of general relativity, given by the action Sg. The Hilbert space 
of all physical states is spanned by the stationary states, with all possible values of E. 

Now multiplying both sides of (50) by NU, where N is an arbitrary function, 
integrating over space, and applying the supermomentum constraint fl93|), we have in 
place of © 

5 2 



1 



J&xJgNU 



d 3 x NG ab (-h 



5q a 5q h 



where 
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+ N i Hl 
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(95) 



5q a 5q b 

is the operator form of the Kinetic/Potential ratio (|29| ) introduced in the previous 
section, with 

N(x) 



N(x) 



mp- 



(96) 



Jd 3 x'^gNU(q) 

It should be remembered that certain operator-ordering contributions to (0), coming 
from the T and T\ terms in (^) have been dropped. However, operator ordering 
terms are always 0(h), and will not affect the correspondence of the JE operator to the 
Kinetic/Potential ratio fl29|) in the classical limit. From the fact that physical states 
are independent of N and Ni 

6 „, 5 



-my 



-JE^f = 



6N(x)~ 5Ni(x) 
we find, as in the previous section, the generalized constraints 
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The closure of these generalized constraints under commutation depends only on the 
assumed closure of the standard constraints (p£2|), for any fixed value of the parameter 
£. 

From the Schrodinger evolution equation (^), we have 



d T < q a > 



— < [JE, q >=< — — > + op.-ordermg terms 

k dp a 



mp 



< I d 3 x' 



T <p a > 



5p a (x) 

+ op. -ordering terms 
-<[JE, Pa ) >=<~ Wa > 



5p a (x) 



> 



op. -ordering terms 



< / d 3 x' 
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5q a (x) 



^G ca p cVd \ 
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+ op. -ordering terms 



> 



(99) 



Equation ( |9~9D is the Ehrenfest principle obtained from our transfer matrix formal- 
ism. Removing the "<>" brackets, replacing JE by a constant parameter —£, and 
dropping operator-ordering terms, eq. ( |9"ED and ( ]9"9| ) become the classical equations of 
motion (|32|) . 



4 WKB, Time, and Many- Fingered Time 

In classical general relativity, the split of spacetime into space + time can be accom- 
plished in infinitely many ways, which are distinguished by a choice of lapse and shift 
functions. The geometrical meaning of a time parameter t, in the classical equations 
of motion 



d t Q = {Q,H} 



d 3 x 



N{x){Q,H x } + N^iQ^} 



(100) 



is specified by N(x) and Ni(x). This is also the case for the iE-evolution equation 



d T Q = {Q,M(p,q,N,Ni)} 
1 



mp 



d 6 x 



y/SN(x){Q,Hi} + Ni(x){Q,K} 



;ioi) 



except that in the latter case the lapse N has been normalized to satisfy the condition 
(|33|). To see what this restriction means for the evolution parameter r, let us consider 



on 



Hamilton's principal function S% P [q, q'\ satisfying the Hamilton- Jacobi equations and 
constraints 



K 2 G a b 



5S £ HP SS £ HP +g ^_ u = Q 



Sq a 5q b 



Pa 
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(102) 



and 



<5</ a 5q" 



Pa 



SS £ HP 



(103) 



As usual, the principal function S HP defined in this way has the interpretation of being 
the action of the classical trajectory in superspace, connecting configurations q a (x) and 



q' a (x) 



Then consider the variation of S HP in time r according to ( |101| ) 



AS' 



dS £ 



HP 



HP 



At = {S £ hp ,M}At 



dr 

d 3 x HP - At 



Sq a Sp a 
5S £ 
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— / (/"./• — — 

m P J oq a 
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$Shp 
5q a 



At 



-2£At 



(104) 



What eq. ( |104|) establishes is that, in evolving classically from configuration q' to 
q in time At according to (|101|) , the action S HP [q, q'\ is proportional to At regardless 
of the choice of N, Ni. In other words, while the particular configuration q reached 
after At does depend on N, Ni, the increment of action AS" does not. This situation 
is indicated schematically in Figure 1. Of course, the restriction At oc AS hp is not 
a restriction on hypersurfaces; Q[q,p] can be computed from ( |101| ) on any spacelike 
hyper surface of the classical manifold, by choosing an appropriate N, Ni. For this 
reason, the iE-evolution equation ( |101| ) is just as informative as the traditional evolution 
equation (|100f) . 

The quantum equation of motion for an observable Q, in our formalism, is 



d T < Q >= -- < V(q,T)]\lQ,m]\*(q,T) > 



(105) 
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The correspondence between ( |105| ) and (|101|) in the semiclassical limit can be studied 
using the WKB approximation. In this approximation, to leading order in h, a solution 
of iE$ 



-£$ and all other constraints (f42|) has the form 



®£,q'(q) = ex P iVSS H p[q,q']/h 



where 



c _ °HP 
JHP = 



(106) 



(107) 



is Hamilton's principal function for £ = 1. Of course, eq. ( |106|) is only valid away from 
caustics, in the classically allowed region. We can now write a general r-dependent 
WKB solution SIS db wavepacket of the form 



*(g,r) 



dSDq' F[S,q'} exp[iS] 

d£Dq' F[£, q'\ exp \i{£t + VSS H p[q, q') + 



(108) 



where F[£, q'} and 9(q') are real- valued. To make contact with classical physics, suppose 
that F is sharply peaked around some value £ and configuration q^x), and define 



M a ( x ), 



(109) 



=90 



Then, according to the principle of constructive interference 0, the wavefunction 
ty(q,r) will be peaked around configurations q a (x) such that the phase S in ( |108| ) 
is stationary w.r.t. small variations of the parameters £ , q' around £ ,qQ(x). In other 
words, q a (x) satifies 




£=£o 



=90 







POa 




(110) 
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The principle of constructive interference is the standard way to make the connec- 
tion between WKB wavefunctions and the classical equations of motion ||, JJ. Equa- 
tions ( P-lOj ), with Shp satisfying the Hamiltonian-Jacobi equations (|102|) and ( |103|) , 
are sufficient to specify a classical solution q a (x, r) beginning from an initial configura- 
tion q a (x, 0) = q$(x), and initial momenta po a (x) satisfying the appropriate constraint 
equations (|103|) . 

The first of the equations in ( |110| ) tells us that, for WKB wavefunctions of the form 
( |108j ), the "wavefunction of the Universe" is not peaked around a single 3-manifold 
at time r, but is instead peaked at configurations (3-manifold + fields) q a (x, r) in 



superspace such that the classical action between the initial configuration q and q a is 
a constant, i.e. 



T 



SHp[q,qo] 



2£ c 

Therefore, even in the classical limit, r does not in general single out a particular 3- 
manifold and set of fields, obeying ( |1 1 1|) . The second equation in ( |110|) is equivalent to 
the dynamical equations of the classical theory, and ensures that the wavefunction is 
peaked only at configurations q a obtained from spacelike slices of the classical solution 
determined by an initial {q^Poa}- 

The situation for F[£, q'} concentrated at some So and go is summarized in Fig. 
2. At a given value of the evolution parameter r, the wavefunction of the Universe is 
peaked at a whole class of configurations, which correspond to spacelike slices of the 
classical, 4-dimensional solution, satisfying the restriction (|111|) . 

Now after a measurement of some observable, it is expected that the Universe 
should be left in an eigenstate, or approximate eigenstate, of that observable immedi- 
ately after the measurement. This is a great difficulty for the standard formulation of 
canonical quantum gravity, since an eigenstate of, e.g., 3-geometry is not a solution of 
the Wheeler-de Witt equation, and therefore not a physical state. In our formulation, 
the situation is more favorable. Imagine that some measurement were performed which 
determined the values of the observables {q a ,p a }, or some subset of these observables 
(modulo three-dimensional diffeomorphisms) up to the limits imposed by the AqAp 
uncertainty principle. Immediately after the measurement, the Universe would be left 
in a state of the form ( |108| ) with r = 0, and go = qobserved (modulo three-dimensional 
diffeomorphisms). This means that, at r = 0, the wavefunction of the Universe would 
be peaked around configurations q a satisfying 

SHp[q, qobserved] = (H2) 

Since: (i) the action SHp[q,q'} of a classical solution bounded by q and q' is only 
well-defined if q' is a Cauchy surface for q (i.e. if q and q' do not intersect on the 
4-dimensional solution manifold); and noting that (ii) the classical action is monotonic 
in the evolution parameter (see eq. (|104| )), away from caustics/turning points; eq. 
( |112j ) implies that the wavefunction at r = is peaked only at configurations which 
are equivalent, up to three-dimensional diffeomorphisms, to 

q a (%) = qobserved(x) (H3) 



In other words, there exist physical states which are at least approximate eigenstates 
of quantities such as 3-geometry, and this is made possible by a dispersion in the value 
of £. 



The "Energy-Time" Uncertainty Relation 



As discussed at length in section 2, the quantity £ is classically irrelevant: it does 
not appear in the Euler-Lagrange equations, and therefore cannot be determined from 
a classical orbit in superspace. In fact, it is easy to see that an uncertainty in £ can 
be interpreted as an uncertainty in the effective value of Planck's constant. This is 
because in the constraint equation (|93) £ can simply be absorbed into a redefinition 
of h, i.e. 

r 6 2 ^ 



ab 



Sq a Sq l 



where 



h 



eff 



■ + ^gUj $ £ = 

h 

7£ 



(114) 



(115) 



Whatever the correct operator-ordering may be, the absorbtion of £ into h is always 
possible. 

We would like to have a quantitative estimate relating the uncertainty in £ (or 
h e ff) to the spread of wavefunctions along a classical trajectory (in superspace) around 
qobserved( x ) > leading to an analog of the time-energy uncertainty relation in ordinary 
quantum mechanics. The discussion below may serve to illustrate some of the issues 
and ambiguities involved in that estimate. 

Consider a WKB wavepacket with F[£, q'\ of the form 



F[£,q] = ^-f[q')exp 



(p - pqY' 

(Apf 



;ii6) 



where p = \/£. Then, integrating over £ in ( |108| ) we have, at r = 0, 



*(<Z,0)= Dq'f[q']e 



o i(e[q']+p S H p[q,q'])/h 



exp[-(ApyS 2 HP [q,q']/4h 



(117) 



Assuming, as before, that f[q'} is peaked (modulo diffeomorphisms) at some configura- 
tion q^ix), the integrand has a stationary phase along the classical manifold statisfying 



the second of equations (|110|) . Consider a configuration q d along the classical manifold. 
Away from q c i = q Q , there will be a suppression factor in the wavefunction 



|^(g d ,0)| 2 ~ exp[-(Ap) 2 ^ P [g d ,g ]/2ft 2 ] 

(^Sgp(<fa,<ft,)) /8h 2 



cxp 
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If we ask only for the spread of the wavefunction along the classical manifold, then 
clearly 

AS 

-^S§ P (q cl ,q Q )~h (119) 

or, taking into account eq. (|104j) , 

AS At ~ h (120) 

which has the form of the usual energy-time uncertainty relation. However, since r is 
just an evolution parameter rather than, e.g., a proper time interval, this relation is 
not very informative. 

Let us consider the probability density P(AS) for being at any configuration q c i 
such that 

AS = S HP [q d ,q } (121) 

In that case, the measure fi(AS) of such configurations along the classical solution 
manifold becomes important: 



P(AS) ~ n(AS) exp 



AS 



SE- 



AS /8h 
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In the absence of a regulator for general relativity, we don't really know the measure 
n(AS). But let us imagine that it has some simple power-law dependence on the 
number of degrees of freedom N p 

fi(AS) ~ (AS) aNp (123) 

where a is a constant of 0(1). In that case, P(AS) would be a maximum at 

/AS \ 2 

AS) = 4ah 2 eff N p (124) 

To relate AS to an increment of proper time As along the classical manifold, let us 
consider a matter-dominated Friedman universe, in which 

AS ~ MAs (125) 

where M is the total mass, and As is the average increment of proper time between go 
and q c i in a synchronous coordinate system. Then we have 

[AS \ 2 - N p 



v _ As j ~aK ffJ £ (126) 
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Let V be the volume and p m the density of matter in the Universe, and denote by v 
the volume per degree of freedom. Then 



AS h eff I a 

—As —J— 77 127 

S Pm V v oV 



or, equivalently, 

Ah , , 

- AS ~ -T-T7 (128) 




which is the analog, in our formulation, of the energy-time uncertainty relation of 
non-parametrized theories. 



It is difficult to assign a reliable number to the right hand side of expression ( 128 ) , 
even assuming the validity of eq. (|123|) as a measure of configurations, which was 
used in the derivation. The mass density p m can be taken as, roughly, the critical 
density for closing the Universe (~ 10~ 26 kg/m 3 ), and Vq, assuming it is different from 
zero, is presumably the Planck volume (~ 10~ 105 m 3 ). But there is still tremendous 
uncertainty in the volume V of the Universe. One could even speculate that, in an 
appropriate continuum limit of regularized quantum gravity, v o —>■ and V — > oo such 
that VoV tends to a finite constant. In any case, while it appears that the transfer 
matrix formulation must generate some uncertainty relation between proper time and 
the fractional uncertainty in Planck's constant, we are so far unable to say anything 
quantitative. 



5 Non-Stationary States in Minisuperspace 

We will now illustrate the formalism developed above in the simplest context possi- 
ble: the quantum mechanical time-evolution of a homogeneous, isotropic Friedmann 
universe filled with a relativistic perfect fluid. In this toy model it is possible to "do 
everything" ; i.e. to find the integration measure and operator ordering, to solve for the 
spectrum and exact eigenstates of the M evolution operator, to form wavepackets with 
a conserved norm, and to track the evolution of such "wavef unctions of the Universe" 
through collapse towards the singularity, "bounce", expansion, and recollapse. 
The Friedmann universe is described by the metric 

ds 2 = a 2 [-N 2 dt 2 + a 2 dtt 2 3 ] (129) 

where N and a are the lapse and scale factors and a 2 = The perfect fluid content 
of this homogeneous and isotropic universe is specified by the energy-momentum tensor 

T^ = {p + p)U^U u +pg^ (130) 



9fi 



where is the Friedmann metric corresponding to the line element (|129|) , = 
(1, 0, 0, 0) with U^U^ = — 1 is the four velocity of the fluid and p and p are, respectively, 
the energy density and pressure of the fluid. 

For a perfect fluid described by the equation of state 



P = (7 - l)p 



:i3r 



the action for gravity and matter (with zero bare cosmological constant) which is 
consistent with the definition of the energy tensor ( |130| ) (see, e.g., [TIJ) can be written 

as 



s> 



bin 



(It 



aa 



+ Na- Na^a 3 ^ 



[132) 



Classically, these models represent, for 7 > 2/3, a universe which expands up to a 
maximum radius a and then recollapses towards the singularity at a = and, for 
7 < 2/3, an inflationary universe expanding from the minimum radius a® (see, e.g., 

H). 4 

The perfect fluid does not introduce any extra dynamical field in this model, and 
it is easy to rewrite the action (|132| ) in the Hamiltonian form corresponding to the 
time-parametrized theory of eq. flT2|) . The only degrees of freedom are given by the 
scale factor (q a = a) and its conjugate momentum p a = — ^7, and the Hamiltonian 7~t 
is given by eq. (|12|) with 



N ' 



Gaa 

V(a) 



2-37 
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(133) 



(and m = 1). The 'superpotential' V is positive definite in the classically allowed 
region a < ao for 7 > 2/3 (a > ao for 7 < 2/3), and it changes sign at the classical 
turning point a = a . 

Following the discussion of section 3, we can immediately write the expression for 
the evolution operator JE in the 1-d minisuperspace (see eq. fl5T|)) as 

h 2 



f = - \g\ 1,2 g aa - 

\Q\ X I 2 da da 



(134) 



In this case there is no operator ordering term involving 1Z since, obviously, a one 
dimensional superspace has vanishing scalar curvature 71 = 0. More explicitly, using 
Gaa — G — — 2aV(a) and eq. (|133[) , we can write 



M = ^ 2 [2a|y|]' 



-1/2 



d_[2a\V\\^_d_ 
da 2aV da 



4 Also note that we must have 7 < 2 for the sound wave velocity of the fluid to be less than the 
speed of light. 
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9a 2 2a y V da J da 
The problem is now to solve for all the stationary states 

m$ e [a] = -£$ e [a] 



(135) 



(136) 



each of which is a solution of a Wheeler-de Witt equation with a particular operator- 
ordering, and an effective value of Planck's constant which depends on £ (see eq. ( jH5p ). 
With these solutions in hand, we can then construct non-stationary states as a linear 
superposition 

^[a,r]=J2^s[a}e l£T/h (137) 
s 

Since the evolution operator M of (|135|) is very singular at the classical turning 
point a = a (the second line of eq. (|135|) does not even make sense, since the product 



of two distributions is not defined) we shall consider the regions where V > and 
V < separately, and then impose some appropriate junction condition at the classical 
turning point a = a . 

Let us first consider the region V > 0. In this region we change from the coordinate 
a to the 'tortoise' coordinate 

1/2 



± da [2Vd} 1/2 = ± dad 

Ja J ao 
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(138) 



where the plus sign is for the case 7 < 2/3 and the minus sign for 7 > 2/3. The 
coordinate a + is chosen such as to start from zero at a = a and to be semi-positive 
definite, monotonically growing to infinity as a > ao — > 00 (7 < 2/3) or to the maximum 
a +M > as a — > (7 > 2/3): 



a + G [0,oo) , a > a , 7 < 2/3 
a + G [0,a +A f] , a < a , 7 > 2/3 
a+ > , a+(a ) = , a+(0) = a +M 



(139) 



With this choice, the eigenvalue problem for the evolution operator simply becomes 



da 1 . 



(140) 



which is the Schrodinger problem for the motion of a free particle with energy £. 
Assuming from now on that £ > 0|] the general solution of eq. ( 140 ) in the region 
V > is clearly a combination of plane waves 



$+[ a+ ] = A e 1 ^^ + B e -^ a +l n 
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3 The choice £ > is required in order to have a real-valued Planck's constant. 
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Similarly, in the region V < 0, we transform to the 'tortoise' coordinate 

. . 1/2 

± / da \-2Va] 1/2 



± / da a 
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(142) 



where, again, the plus sign is for the case 7 < 2/3 and the minus sign for 7 > 2/3. 
The coordinate a_ is zero at a = ao and negative semi-definite, monotonically growing 
to minus infinity as a > a — > 00 (7 > 2/3) or to the minimum a_ m < as a — > 
(7 < 2/3): 



a_ G [a_ m , 0] , a < a , 7 < 2/3 

a_ G (—oo,0] , a > a , 7 > 2/3 
a_ < , a_(a ) = , a_(0) = a_ 

In this case, eq. ( |136| ) becomes 



q2 



with solutions in the (classically forbidden) a_ region 

$j[a_] = C e ^ a - /n + D e -^ a - /h 



(143) 



(144) 
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Let us now build the general solution of the original eigenvalue problem (|136[) in 
the whole range a G [0, 00). 

First, we must impose a junction condition at the turning point a + = a_ = 
[a = ao J. We multiply both sides of eq. flUKf) by [2a|V|] 1/2 and, using eq. (|135|) , we 
integrate around the turning point ao- Assuming a reasonably smooth behaviour of 
at ao, we find that 

f a o+<5 . . ..„,„£ 



ro+s \2a\V\] l l 2 d 
lim / da — — — — 

s^oJa -8 oa 2aV da 



$ £ = - lim / da [2a|1/|] 1/2 -2$ f = (146) 



Evaluating the integral on the left hand side, the joining condition on first derivatives 
of $£ in the a coordinate can be rewritten as 



Wl 1 / 2 d$ f 



V da 



l 0+ 



V da 
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Both $£■ and its derivative must be continuous at a = ao. 

When we turn to a± coordinates, it is easy to see that first derivatives must be 
'discontinuous', i.e. we have the junction conditions 



*+[0] = *J[0] 



daj 



da^r 



[148) 
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where in the last equation upper signs refer to the case 7 < 2/3, and lower signs to the 
case 7 > 2/3. 

Second, the evolution operator is defined on the 'half line' a > 0, and therefore it is 
not essentially self-adjoint in L 2 [0, 00) with the measure \f\G~\- One can build, however, 
a one-parameter family of self-adjoint extensions (which guarantee norm conservation 
and unitarity) by appropriately choosing boundary conditions at the origin (see, e.g., 
flUl). In the a± coordinates this is translated into the condition 



da± 



a< &£ \ a±, m (M) 



(149) 



a ±,m(Af) 

where a is an arbitrary real parameter with range (—00, 00). The lower signs (and index 
m) correspond to the case 7 < 2/3, while the upper signs (and index M) correspond 
to 7 > 2/3. 

Self-adjointness of the evolution operator M automatically guarantees orthogonality 
of eigenfunctions, i.e. 

c0 



da \j\G\§* s A\a\<&s 2 V 



+ 



a_ m (— 00) 
oo(a +M ) 



da- $ £l [a_]$ & [a_ 







da + *t*[a + ]$l[a+]=6(£ 1 ,£ t ) (150) 

where 6(81,82) is the Dirac delta function when the spectrum is continuous and the 
Kronecker delta for a discrete spectrum (and integration limits in parenthesis are for 
the case 7 > 2/3). 

It is also possible to show that the system of eigenfunctions is complete, i.e., defining 
p=\J~8, that 

6 (a — a') 
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(where of course the sum is meant as an integral for a continuous spectrum). 

We are now in the position to write the exact eigenfunctions for different 7 as a 
function of the scale factor a. 



• 7 < 2/3 

In the case 7 < 2/3, the spectrum is continuous and non degenerate, and from 
conditions ( |148| )-( |i~4"9D one finds, for a generic boundary condition at the origin a = 0, 



.4 
B 
D 



i-i)^M e ^- 

(V8 + ha) 



(V8 - ha) c2 ^ a 
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For instance, in the case 7 = 0, which classically corresponds to a de Sitter universe 
expanding from the minimum radius a , the 'tortoise' coordinates derived from eqs. 
( Hp and (HD are 



a J 



1 3/2 



a+ 



- 1 



3/2 



Moreover, if we choose boundary conditions at the origin such that $^(0) 
equivalently, from eq. ( |149| ), such that a = oo), the eigenfunctions of M are 



(153) 
(or, 



®e[a} = Cle 
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for a < a and 



2Ce~ 3R <^ sinh 
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for a > a . With this choice of boundary conditions, the eigenfunction is real for all 
values of a. 



• 7 > 2/3 

The case 7 > 2/3 corresponds classically to a universe which expands from the 
singularity at a = to the maximum radius at a = ao, and then collapses back towards 
the singularity. Since now the range of the 'tortoise' coordinate a_ is infinite (see eq. 
( |143| )), to ensure square integrability of the eigenfunctions (eq. ( |150| )) clearly we must 
have D = 0. Therefore, due to conditions ( |14£| ) and ( |149| ) , the spectrum is discrete 
and non degenerate. 

Normalizing the eigenfunctions to one in the range a e [0, 00), from eqs. (|148|) and 
( |150j ) we have that 

a = (i+ik 



31 
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C 



f l-i) 



c 



2V£ 



h cos(2\/£ a + M/h) + 2\/£a + M 



1/2 



(156) 



where 9 is an arbitrary phase. 

Moreover, imposing a generic boundary condition at the origin according to eq. 
( |149| ) gives an implicit expression for the discrete spectrum of £, i.e. 



tan(VSa +M /h) 



ha + \[E 
ha — \[E 



(157) 



For example, in the case of a pressureless, dust-dominated universe with 7 = 1, the 
explicit form of the 'tortoise' coordinates is 
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(158) 



Choosing boundary conditions such that $£-(0) = (or a = 00), from eq. ( |157|) one 
can explicitly write the discrete eigenvalues for the dust-dominated universe as 



v?^ 2 ( 4 "t 1)fi 



n = 0,1,... 



(159) 



and the JE eigenf unctions are 
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for a < ao and 
$e[a] 
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for a > a . 

In the same way, one can solve for the stationary states of a radiation-dominated 
universe with 7 = 4/3. In this case the discrete eigenvalues are 



(An + l)h 



«o 
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and the /E-eigenfunctions are 
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for a < a and 



2e l 
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exp 
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for a > ao. The main properties of the wave function are essentially the same as those 
of the dust-dominated universe. 

Note that in both the 7 > 2/3 and 7 < 2/3 cases, the M eigenstates are real- valued. 
In the more conventional "Born-Oppenheimer" interpretation [16| of the Wheeler-de 
Witt equation, such wavefunctions represent a superposition of 'collapsing' and 'ex- 
panding' universes in the classically allowed region. In our formulation, by contrast, 
such states are stationary, and dynamics (such as the expansion and collapse of the 
universe) arises from their superposition. We will now illustrate such dynamics by 
constructing a non- stationary state from the M eigenstates. 

Let us consider the case of a 7 > 2/3 Friedman universe, which has a discrete 
spectrum for S, and oscillating eigenfunctions in the classically allowed range a G [0, ao]. 
The general form of the wave packet for an arbitrary distribution f[p n ] of "momenta" 
Pn=V£n is given by the standard formula 



a\e 



n=0 



where 



f[ P n}= J Q da 0]$; 
and, for the choice of boundary conditions &s(0) = 0, 

7r(4n + l)h 
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We build up the wave packet such that it is initially localized in the 'classically allowed' 
region a £ [0, ao] and it moves towards the singularity with "momentum" p n = \[E^ 
centered around p no . For our example, we choose an initial wavepacket which is a 
triangle of unit height, with base of length 25 centered at a + o, i.e. 



tf[a+,0] =e 



1 - 



a+o| 



[6(a+ - a +0 + 5) - 6(a + - a +0 - 5)} (168) 



where 9 is the Heaviside step function. 

Since the initial packet has support only in the region a + , from eqs. (|141|) , (|166| ) 
and ( |168|) we find 

r a +M 
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(Pn + Pn) 

Inserting this result in eq. ( |165| ) and taking the modulus square, we finally find that 
the probability distribution of the wave packet evolves, in the classically allowed region 

cz_i_ , as 
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where, for convenience of notation, we have defined 
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Similarly, the probability distribution in the classically forbidden region a_ is easily 
found to be 
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where 
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For our illustrations, we have chosen to compute the probability distribution for 
a dust-dominated universe (7 = 1) with n Q = 400, a = 1 and 5 = 0.02 (keeping 
10 4 terms in the series). The results are plotted in fig. 3, where different times are 
labeled by r*== — 10 5 r. For a comparison, we have also plotted in fig. 4 the modulus 
of the eigenfunction of M for the S eigenvalue corresponding to the average £ of the 
packet, namely £ = (3202) 2 . This stationary state is the solution of a Wheeler-de 
Witt equation with a particular choice of operator-ordering and Planck's constant. 

Fig. 3 shows the evolution of the wave packet; the packet begins by moving 
smoothly towards the singularity while gradually spreading. As it approaches the 
classical singularity at a = 0, the packet starts oscillating, with the frequency and am- 
plitude of oscillations increasing up to a maximum which is shown in more detail in fig. 
5a. The wave packet is always zero at a = 0, consistent with our choice of boundary 
conditions. As r increases, the packet 'bounces' off a = and moves back towards the 
classical turning point at a = a^. At the point of maximum expansion the wave packet 
has again an oscillating behaviour similar to that close to the classical singularity. Here, 
however, part of the wave packet extends into the 'classically forbidden' region a > a ; 
there is a small, exponentially decaying probability to find the universe in such a region 
(see fig. 5b). At a = ao both the wave packet and its first derivative are continuous. 
As the time parameter continues to increase, the wave packet 'bounces' off a = ao and 
moves back towards the classical singularity. The packet rebounds repeatedly between 
a = and a = a$ while gradually becoming delocalized. It is important to stress (see 
also ||) that nothing pathological happens to the wave packet as it bounces between 
the classical singularity at a = and the turning point at a = ao. In particular, there 
is no infinite compression of the wavefunction, and the transition from 'big crunch' to 
'big bang' is a (relatively) smooth process, in this toy model, at the quantum level. 

The norm of the wave packet has been numerically checked to be constant during 
all the stages of the evolution, as it obviously should. 

Although we have specialized here to the case of a dust dominated universe, it is 
quite clear that we do not expect anything radically different from the analysis of other 
cases with 7 > 2/3. The same holds for the cases 7 < 2/3, the only difference being 
that there is only a single bounce of the kind shown in fig. 6, with the wave packet 
starting from infinite scale factor, 'bouncing' off the minimum radius ao and finally 



going back to infinite scale factor. Also, of course, the qualitative picture does not 
depend critically on the shape of this initial wavepacket. 

As pointed out in section 2, the space of physical states is spanned by eigenstates of 
the evolution operator M with different £, and in section 4 it was shown that this can 
be interpreted as leading to a quantum indeterminacy of Planck's constant according 
to formula ( |115| ). In our minisuperspace toy model, it is possible to directly relate the 
spread of the wavepacket in the scale factor a to the spread in S, and thereby to the 
dispersion in Planck's constant. Of course, this model is far too unrealistic to draw 
any quantitative conclusions regarding the dispersion of fundamental constants vs. the 
spread of the wavepacket in full quantum gravity. 

To proceed, we construct a wave packet from eigenf unctions of M. Now, however, 
we want to introduce a more tractable Gaussian distribution in the momenta, and to 
simplify the analysis let us assume that the £ eigenvalues are so closely spaced that we 
can replace sums over the discrete eigenvalues by an integral,^ i.e. we write 

dpf[p]$ p (ay^ h 
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-P 2 (P-P0) 2 
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where 



vjHWr (177) 

is the dispersion in the "momentum" p = \f£ distribution around the peak at po- 

Evaluating the integral in the classically allowed region with $ p = e~ ~ and taking 
the modulus, we easily find that 



\fy[a, t}\ 2 ~ exp 
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2 h 2 
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where a + (a) is given by formula (|138|) . 

Following standard analysis (see, e.g., |T7[), we see that, at any given r, the probabil- 
ity distribution of the scale factor a is peaked around the set of semiclassical trajectories 
for which 

a + {a)\ d = -2sf£~ t (179) 

or, in other words, for which the phase in the integral ( |176l) is stationary (compare 
with eqs. ( TO and (|TT0D). 



As an example, we can differentiate both sides of eq. ( |179|) with respect to r and 



find 
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3 This can be done by replacing X^^Lo ™* IT^F I^oo f° r 7 > 2/3. 



If we now compare with the known solutions in classical cosmology for the perfect fluid 
matter content (see |T5[), i.e. 
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we can immediately infer the correspondence between the 'internal' time t and the 
evolution parameter r at the classical level, given by 
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o_\ 2 -^ _ 
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We now expand |\l/[a,r]| 2 in eq. ( |178|) around the classical trajectories with a 
a(a + \ c i) + Aa at r fixed, and find 



( a ~ a cl> 



where 



T 



h 2 (3 2 



da 



2 (2hG\ 2 



3tt 



and we have reinserted dimensions according to 
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Eq. 



can be interpreted as saying that the probability to find the universe in a 



certain configuration with scale factor a at time r is given by a Gaussian peaked at 
the classical trajectory a d with a dispersion a 2 a . 

Finally, as an exercise, we can tentatively give a lower bound on the dispersion for 
the scale factor by minimizing it with respect to (3, and we find 



{Min[a 2 a }} MT/nl i/2 = |a H 



da + 
da 



T 2 



(186) 



and the 'dressed' Planck's constant h e ff is 



where we have defined L 2 Pe ^ - .,_ 
introduced according to eq. ( |115|) . We can study the behaviour of eq. ( |186| ) for some 
simple examples of matter content in the universe. In the case of a de Sitter universe 
with 7 = we can use formula ( p.53[ ) for a + , define a = £a (£ > 1) and find 



{Mtn[a 2 a }} MT/nll/2 _ (e - iy/ 2 ^ 1 
Lp, e ff 3£ 2 36 



V^> 1 
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Therefore, if the dispersion in the scale factor is such to be minimized with respect to 
(3, it will always be negligible for the case 7 = 0. A similar conclusion can be easily 



^7 



drawn by studying the case of a radiation-dominated universe. Unfortunately, this 
simple analysis does not either provide any reliable upper bound on a a . 

For a dust-dominated universe, instead, using eq. (|158|) with a = (,a (£ < 1), one 

has 

[arcco 3 (0'/ 2 + (l-20K(l-0]'/'] 

In this toy model, it is therefore possible to have non-negligible dispersion close to the 
classical singularity, while having a a negligible at much later times. 



6 Conclusions 

In this article we have pointed out that there is a constant (denoted £ ) in the classical 
Hamiltonian of certain parametrized theories, which cannot be determined from the 
classical trajectories. It is the arbitrariness of this constant in classical physics which 
allows us to extend the space of physical states in the corresponding quantum theory, 
and this is the key to our proposed resolution of the problem of time. Solutions of 
the Wheeler-de Witt equations are stationary states, in our scheme, but there is an 
infinite set of such equations, each with its own parameter £. Non-stationary states 
are composed of superpositions of such states, with different £ parameters. It was 
seen that this formalism can be derived directly from a transfer-matrix quantization of 
parametrized theories. 

As noted in section 4, the constant £ in the Wheeler-de Witt equation (equivalent to 
= —£&£ in our notation) can always be absorbed into a redefinition of Planck's 
constant h e ff = Uj\f£. In the special case of pure gravity, the constant £ could 
alternatively be absorbed into a redefinition of Newton's constant ( (4"5| ) (as was done in 
ref. j|]), but in pure gravity this is anyway equivalent to rescaling the Planck constant, 
since only the combination hG/\/£ appears in the Wheeler-de Witt equation in this 
case. In the general case, any non-stationary state 

nq,r]=J2^s[q]e i£T/h (189) 

£ 

involves a superposition of Wheeler-de Witt wavefunctions $£•[(?], each with a different 
effective value of Planck's constant, so in this sense any non- stationary physical state 
entails an inherent uncertainty in the value of Planck's constant. Unfortunately, we 
have not been able to place a reliable lower bound on this uncertainty, for reasons 
which were discussed in section 4, above. Nevertheless, the possibility that there is an 
inherent uncertainty in the effective value of Planck's constant, which might be large 
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enough to be observable, raises some interesting phenomenological questions, which we 
hope to address in the near future. 
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A Appendix 



We briefly comment here on a 2-d minisuperspace model for a universe filled with a 
homogeneous scalar field <fi with arbitrary potential V(0) and with a bare cosmological 
constant A. The geometrical ansatz is the same as that of the 1-d minisuperspace 



model, i.e. the isotropic and homogeneous Friedmann metric given by eq. (|129|) . The 
reparametrization invariant action for such a model is 
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This can be easily written in the Hamiltonian form of eq. ([)!]), with p a 
m — 1 and 
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The general form of the evolution operator JE is 
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where, as usual, we have introduced the 'supermetric' Q mn = 2VG mn , with determinant 
Q and corresponding scalar curvature 
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As a consequence of the hermiticity requirements for the measure, similarly to the 
1-d minisuperspace case, the formal expression of the evolution operator is not well 
defined at the 'turning points' where V = 0. In fact, from eqs. (|191|) -( |T93| ) and using 
the explicit form of the supermetric, one gets 
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The general procedure will be to study the eigenvalue problem for M (eq. (|136|) ) in 
the two regions where V > and V < separately, and then impose suitable joining 
conditions on the partial derivatives of the eigenf unctions at the 'turning point' V = 0. 

The simplest problem which we can analyze is that corresponding to a massless and 
minimally coupled scalar fieldj] with 



V(cf>) = A = 



(195) 



for which V — — | and TZ — 0. In this case, since the scale factor a is constrained to 
be positive, the evolution operator given by eq. (|194|) has no ambiguities due to delta 
terms and becomes 
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The hyperbolic problem described by eq. ( |136| ) with the operator ( |196| ) and with 
£ = 1 has been extensively discussed in the literature as an ordinary differential prob- 
lem (see, e.g., [|18[ and references therein) or as a Klein-Gordon eigenvalue problem 
with respect to a peculiar measure (see, e.g., [0). One can separate the differential 
equation and express bounded solutions as a superposition of modified Bessel functions 



of imaginary order |22 



e,k[ a i 



exp (±ik<p/h) Kik_ 



VSa 2 \ 



(197) 



where k is the separation constant, or give them as a product of harmonic oscillator 
wave functions like 
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where X = a cosh 6 and Y = a sinh 6 ITS 



Unfortunately, neither of these solutions are orthogonal with respect to the measure 
\J~G oc a 3 , as it is easy to check. To guarantee the self-adjointness of the evolution 
7 Thc case of a conformally coupled scalar field turns out to be essentially equivalent. 



in 



operator, one has to construct a particular superposition of any of such solutions which 
satisfies appropriate conditions on the boundary of the 2-d minisuperspace. 

A possible method of solution of the eigenvalue problem makes use of the so called 
'Rindler' coordinates pl| . Here we derive the formal boundary conditions for the 



eigenfunctions required for Hermiticity, even though we have not been able to explicitly 
solve for them. The idea is to make the following transformation of coordinates in the 
2-d minisuperspace : 

x = cosh 26 
V6 

a 2 

y = -^sinh20 (199) 
V6 

The minisuperspace line element in the new coordinates takes the simple form 
ds 2 Super = (dy 2 — dx 2 )/4. The minisuperspace region a 6 [0, oo) and 6 e (—00,00) is 
mapped into the right Rindler wedge of a 2-d Minkowski spacetime, with boundary 
\x\ < y. The correspondence between old and new coordinates is one-to-one, with 
a = mapped to x — ±y, 6 = constant to - = constant lines and a = constant 
mapped to the hyperbolae x 2 — y 2 = constant. The evolution operator simply becomes 
the 2-d wave operator in the Rindler wedge, i.e. 

R - fi2 (w - & (200) 

By further 'tilting' the minisuperspace axes by 45° with the coordinates 

* + - i<* + »> 

x- = \{x-y) (201) 
(x± G [0, 00)), the evolution operator becomes 

M = ~ h2 w^— ( 202 ) 

OX+OX- 

with constant measure. The general solution of the differential equation (|136|) is thus 
a linear combination of plane waves 

= JdX [A(\)e i{Xx+ -^ x ~ )/h + B(\)S Xx —^ x+)/n ] (203) 

Moreover, the self-adjointness condition for the evolution operator M requires that 

rod rod Q% rco rod 

\ d.r . 
Jo Jo 



dx + I dx ~ ^a^SZ ^ = /„ dx + h dx - * & 5ZaZ v « (204) 



1 1 



Integrating by parts, one can see that this is true if the following boundary condition 
for the eigenfunctions holds 



dx- 



a 



dx- 



£2 



dxj. 



* 9 , 



Si 



(205) 



Alternatively, one can transform the 'infinite' boundary to a finite one by use of 
the following conformal transformation on the superspace metric ds 2 Super : 



d? 2 

u " 5 Super- 



it 



il 2 ds 2 Super = —dudv 
2 arctan[(x — y)] 
2 arctan[(x + y)] 



u E [0, 7r] 
u G [0, 7r] 



In this case the evolution operator transforms as 



JE=-16/iW(?) cos*^ 



5 2 



dudv 



with constant measure, and it is self-adjoint if 



J u=0 



(206) 



(207) 



(208) 



w=0 



The eigenfunctions are the same as those given in eq. ( [2031 ), provided one makes the 
substitutions x + — > 2tan|, x_ — > 2tan| and £ — > £/16. 

Thus, formally, the hyperbolic eigenvalue problem of eq. ( |136|) for the 2-d min- 
isuperspace model with the minimally coupled scalar is solved provided one fixes the 
coefficients A and B in the expansion ( |203| ) such that to satisfy given Cauchy data and 
the condition ( |205| ) (or ( |208|) ) at the boundary In particular, although we cannot give 
an explicit formula, the problem ( |208j ) certainly admits a solution. A possible simple 
example would be, for instance, a wave function such that $>s(u, 0) (3>g(0, v)) has a 
fixed parity with respect to the point u = § (v — §). 
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Figure Captions 



Fig. 1 Classical evolution via d T Q = {Q, M}. The final configuration q', q", ... depends 
on the N, iV, functions, but the classical action S £ [qj, q ] from the initial to final 
configurations is proportional to the evolution parameter r, independent of N, iVj. 

Fig. 2 Expansion of the (semiclassical) wavepacket. WKB wavefunctions ty[q, r] arc 
peaked, at any parameter time r, on a set of configurations {q c {\ corresponding 
to hypersurfaces of the classical 4-manifold, such that r = —S £ [qo,q c i\/2£- 

Fig. 3 Evolution of a 'triangular' wavepacket for a dust dominated universe. The 
parameters chosen are n = 400, a — 1, 8 — 0.02 and different times are labeled 
by n = -10 5 r. 

Fig. 4 The modulus of the eigenfunction of M for the average eigenvalue So = (3202) 2 
of the wave packet. 

Fig. 5 a) The wave packet at its closest approach to the classical singularity, 
b) The wave packet at the radius of maximum classical expansion. 
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time T ^0: ¥[q,T ] 
is peaked at q=q 
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